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We present a four-body semileptonic charm decay D + — > K~ir + fi + v analysis in the range of 0.65 GeV/c 2 
< rn-Kir < 1-5 GeV/c 2 . We observe a low mass scalar contribution of 5.30 ± 0.74 lo'gg% with respect to the total 
D + — » K~ 7r + jj, + v decay, compatible with the phase shift found by the LASS elastic scattering experiment. For 
the K*(892)° resonance, we obtain a mass of 895.41 ± 0.32 ± MeV/c 2 , a width of 47.79 ± 0.86 t ^ MeV/c 2 , 
and a Blatt-Weisskopf damping factor parameter of 3.96 ± 0.54 o go GeV -1 . We also report 90% CL upper limits 

of 4% and 0.64% for the branching ratios F ^(^^K- 6 ^T^+t)^ and ^wgT^^=lf^7r^T^ ' res P ec tively. 



1. Introduction 



Weak semileptonic decays of charm mesons 
continue to attract interest due to the relative 
simplicity of their theoretical description: the ma- 
trix element of these decays can be factorized as 
the product of the leptonic and hadronic currents. 
This makes the D+—>K~ "7T + n + v decay a natural 
place to study the Ktt system in the absence of 
interactions with other hadrons. 

It is known that the Ktt final state of D + — ► 
K~tt + [i^v decay is strongly dominated by the 
K* (892)° vector resonance |1I2| . The large and 
clean sample of D + — >K~tt + n + v events collected 
by the Fermilab FOCUS experiment provides an 
excellent opportunity to measure the K* (892)° 
mass and width, as well as the effective Blatt- 
Weisskopf damping factor parameter discussed in 
Rcf. |31. 

In addition we also search for structures other 
than the if* (892)° resonance in the mass range 
of 0.65 GeV/c 2 < m K7r < 1.5 GeV/c 2 . The 
available phase space allows us to investigate the 
low mass scalar sector whose spectroscopy has 
been increasingly studied since the indication of 
the existence of the broad scalar resonance k in 
D + —>K ~tt + tt + decays from a Dalitz Plot analy- 
sis performed by E791 g]. Recently E791 has also 
pursued an independent approach for the scalar 
component extracting directly the scalar phase 
shift from their data pj- These two results, ob- 
tained using D + — > K ~tt + tt + hadronic decay is 
substantially different from the one obtained in 
a partial wave analysis performed by the LASS 
scattering experiment 0. LASS observed that 
the s-wave amplitude in the K~tt + elastic region 
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Figure 1. Comparisons between the results ob- 
tained for the s-wave Ktt scattering phase shift 
measured by LASS and the one measured by the 
E791 experiment (solid circles) in D + — >K~tt + tt + 
charm hadronic decay (Figure provided by B. 
Meadows) . 



can be represented as the sum of a ^(1430)° 
resonance and a smooth shape, consistent with 
the non-resonant hypothesis. Figure ^ shows a 
comparison of the phase shift of the scalar com- 
ponent found by LASS to the one found by E791 
in D + -^K~it + tt + decay. 

Due to Watson's Theorem |7l8j the observed 
Ktt phase shifts in D + — > K~tt + h + v decays 
should be the same as those measured in Ktt elas- 
tic scattering. This establishes a connection be- 
tween the production and the scattering mech- 
anisms, adding crucial information about Ktt 
hadronic interaction and the role of the bachelor 
tt final state interaction in hadronic decavs. 
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2. The Experiment and the Data Sample 

The data were collected in the Wideband pho- 
toproduction experiment FOCUS during the Fer- 
milab 1996-1997 fixed-target run. In FOCUS, 
a forward multi-particle spectrometer is used to 
measure the interactions of high energy photons 
on a segmented BeO target. The FOCUS de- 
tector is a large aperture, fixed-target spectrom- 
eter with excellent vertexing and particle iden- 
tification. The FOCUS beamline 0, the detec- 
tor |10I11I12I13| and the specific cut selection ^1] 
have been described elsewhere. A total of 18245 
D + K~n + ii + v candidates remain after the se- 
lection criteria. We estimate the charm and non- 
charm background contributions to be, respec- 
tively, 17.8% and 3.2% of the total number of 
events over our signal region. 

3. The D + — > K~ n + signal description 

Four-body decays of spinless particles are de- 
scribed by five kinematic variables. The variables 
chosen in this analysis are the K~n + invariant 
mass (tokw); the square of the ^i + v mass (q 2 ), and 
three decay angles: the angle between the tt + and 
the D + direction in the K~tt + rest frame (6 V ), 
which defines one decay plane, the angle between 
the v and the D + direction in the y^v rest frame 
(6t), which defines the second decay plane, and 
the acoplanarity angle (x) between these two de- 
cay planes. The differential decay rate can be rep- 
resented by a coherent sum of resonant and non- 
resonant contributions to the angular momentum 
eigenstates of the K ~tt + system, 



r/r 



dm 
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where dtt = dq 2 dcos8 v dcosdg dx, M-j is the 
weak matrix element for a transition with angular 
momentum J, Ajjr) represents the form of the 
hadronic final state amplitude contribution of res- 
onance R (or non-resonant) with strength aj^, 
and <j> is the phase space density. 

The possible resonant states that couple to 
K~ir + are the scalars k and ^(1430)°, the vec- 
tors K*{892)° and K*(1680)°. 2 The non- 
resonant contribution is assumed to be scalar. 



2 Due to the orthogonality of states with different angu- 
lar momentum, only amplitudes with the same spin will 
produce significant interference contributions to the m^-n 
mass spectrum. The inclusion of a small i\^(1430)° res- 
onance contribution is unlikely to be observed, since it 
is orthogonal to the (dominant) ftf*(892)° and low mass 
s-wave amplitudes. 



The parametrization of resonant states with 
angular momentum J is given by the product of 
a Breit-Wigner and the normalized R — > K ~tt + 
coupling, Tj 



Aj.R = 



too r 



m K7r -m^ + imo T(m K7r ) 



(2) 

where T(m Kv ) = T T 2 g p* is the mag- 
nitude of the kaon momentum in the resonance 
rest frame, p$ = p*(m ), T§ = 1, and T\ = 
p* " ^ 1S ^ e Blatt-Weisskopf damping fac- 



tor given by B = 1 / y/l + r 2 , p* 2 j2] . The damp- 
ing factor adds an additional fit parameter, ro, 
in our fits to the K* (892)° line shape. The line 
shape of the k resonance is expected to devi- 
ate significantly from a pure Breit-Wigner, due 
to its large width and the close vicinity of the 
Kn threshold. In this analysis we use the n line- 
shape adopted by E791 0] while for the non- 
resonant component we use an empirical phase 
shift, parametrized by the Effective Range model, 
obtained by LASS (Fig. d 3 . 

The weak matrix element for the vector pro- 
cess, Mi, and for the scalar process, A4q, are 
written as a function of helicity amplitudes, Hi, 
derived in |15|. 

4. Angular Distribution Results 

Next we discuss the angular distribution de- 
scribed by Eq.n The Kir spectrum described by 
this equation includes the dominant contribution 
from the K* (892)° resonance, possible high mass 
contributions from the K$(U30)° and K*(1680)° 
resonances, and low mass scalar components com- 
prised of a non-resonant and a possible k contri- 
butions, both populating the region where signif- 
icant discrepancies were found between the data 
and the predicted D+ -> 7?*(892)°/z+i/ angular 
decay distributions. It has been shown in |l()llfi| 
that a nearly constant amplitude and phase con- 
tribution to the helicity zero amplitude of the 
virtual W + was required to adequately fit the 
observed decay angular distributions. The tok^ 
distribution weighted by cos 8 V provides informa- 
tion on the phase of the additional structure rel- 
ative to that of the K * (892)°. It can be used to 
discriminate different combinations of low mass 
states, given the large difference between their 
phase shifts. Figure compares the distribution 



3 We have to remove the two-body phase space factor, 
given by , from LASS non-resonant amplitude, which 

is already included in Ea.ITl 
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Figure 2. The background subtracted distri- 
bution of mK-rr weighted by cos# v - The data 
(squares) show good agreement with the LASS 
non-resonant parameterization (solid histogram) 
but not with a k model (dashed histogram) 

obtained in the data with the predictions from 
the non-resonant and k models in the absence of 
additional phase shifts as required by Watson's 
Theorem. Since a simulation using the LASS 
parametrization of the non-resonant contribution 
is sufficient to reproduce the data, we exclude 
a possible k contribution from further consider- 
ation. 

5. Kir Mass Spectrum Analysis 

Having excluded the K, the most general differ- 
ential decay rate for D + — > K~n + n + v in is 
given by Eq. |3| 



dT 

dm K% 



= |V| 2 F n + |5| 2 F 00 



2 5i(V* S)F 01 (3) 



The Fjf = J eM*jMf 4>dVl, are real func- 
tions that depend only on m^- These functions 
are computed from the toktt spectrum obtained 
from a complete simulation of D + — > K~tt + [i + v 
events, generated according to phase space and 
weighted by M*jM f and thus represent the in- 
tensity modified by acceptance and efficiency, e. 
The vector, V, and scalar, S, amplitudes are given 
by the sum of the relevant states, each one rep- 
resented by the product of a real magnitude, a, 
and A. The |V| 2 , |S| 2 , and K(V*S) functions de- 
pend on ttt-kw as well as on all fit parameters. The 
cross-term, 2 5i(V* S)F i, represents the interfer- 
ence between the vector and scalar contributions. 

The fit parameters are the magnitudes of 
each amplitude, ai, the mass and width of the 

4 A11 imaginary pieces of A4jM ji will appear as sinusoidal 
functions of \- Hence any imaginary terms vanish when 
averaged over \ given our nearly uniform acceptance in 
this variable. 



if* (892)°, and the parameter r of the Blatt- 
Weisskopf damping factor. The parameters of all 
other resonances are fixed to the PDG values |17) . 
The fit parameters are obtained from an unbinncd 
maximum-likelihood method, by minimizing 
the quantity -2 In £ eventB [(1 - f B )C s + JbCb], 
where Cs is the probability density for the sig- 
nal, Cb, for the background, and fs is the back- 
ground fraction, fixed in the estimated values. 
The K*(892)° amplitude is taken as the reference 
amplitude. Decay fractions are obtained inte- 
grating each individual amplitude over the phase 
space and dividing by the integral over the phase 
space of the overall amplitude. 

To account for momentum resolution effects on 
the K* (892)° parameters, we refit the data fixing 
all parameters except the K* (892)° width and use 
a probability density function convoluted with a 
Gaussian distribution with a = 5.88 McV/c 2 , 
value obtained from Monte Carlo simulation. The 
smearing due to momentum resolution increases 
the K *(892)° width by approximately 2 MeV/c 2 . 



6. Kir Mass Spectrum Results 

Using the procedure described above, we fit the 
data assuming only a D + — > K * (892)° is pro- 
cess. The confidence level of this fit is 0.21%, 
indicating the need for additional contributions 
in the decay. 

The inclusion of a non-resonant scalar compo- 
nent, referred to as the NR model, significantly 
improves the confidence level of the fit to 66%. 
We find m K , (m2)0 = 895.41 ± 0.32 MeV/c 2 , 
Lk« (892)0 = 47.79 ± 0.86 MeV/c 2 , r = 3.96 ± 
0.54 GeV" 1 and a NR = 0.327 ± 0.024, which 
correspond to a scalar fraction of 5.30 ± 0.74 
%. Figure 01 illustrates the contribution of both 
the D + — > K* (892)° v and non-resonant s-wave 
process to the observed mKw spectrum. 

We also consider D+ -> K* (1680)° fi+ v 
and D + — > K$(\A?>Q) [i + v contributions to 
our model. The observed _ftT*(1680)° and 
iQi(1430) amplitudes are consistent with zero 
and we find ^^^f^ < 4.0% and 

r(D+^7fg(1430)V+i/) 



r(D+^>K-Tr+v+v) 



< 0.64% at 90% CL. The 



upper limits are calculated using the method 
described in Ql| and assume BR(^*(1680)° -> 
K-7T+) = 0.258 and BR(K*(U30)° ^ K-n+) = 
0.62 17 . To study the statistical significance 
of these new amplitudes, we use a hypothe- 
sis test based on the maximum-likelihood ratio 
method [HJ\. As a result, we obtain a confidence 
level of 80% in favor of the simole NR model. 
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Figure 3. Fit to the data using the NR 

model. The error bars, the solid lines, the dashed 
lines, and the dotted lines correspond to the data, 
the model, the background contribution, and the 
scalar contribution, respectively. The upper right 
plot shows the same information and the cross- 
term (dot-dash line) with a limited y-axis to allow 
more detail to be seen. 

7. Summary and Discussion 

In conclusion we find that our angular distribu- 
tion is consistent with the effective-range scalar 
non-resonant phase shift obtained by LASS 6 
as expected by Watson's Theorem given the ab- 
sence of other final state interactions (FSI). It 
points out that the FSI of the bachelor pion of 
the D — + (Ktt)tt decay can play an important 
role towards the observed difference with respect 
to Kir scalar scattering phase shift. 

We have measured the i^*(892)° parameters 
using a large sample of D + — ► K* (892)° fj, + v sig- 
nal events over a wide mass range. The absence of 
high mass resonances as well as the small back- 
ground contribution provides a unique environ- 
ment to study the K* (892)° mass and width. Our 
measurements of the mass and width are more 
than 1 a below the present world average value. 
We obtain a Blatt-Weisskopf parameter consis- 
tent with the value obtained by LASS [Hj • We also 
limit possible additional Kir resonances present 
in D + — > K~tt + n + v semileptonic decays. 
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